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Abstract 

We present a detailed analysis of the recently demonstrated technique to generate 
quasi-stationary pulses of light [M. Bajcsy et ai, Nature (London) 426, 638 (2003)] 
based on electromagnetically induced transparency. We show that the use of counter- 
propagating control fields to retrieve a light pulse, previously stored in a collective 
atomic Raman excitation, leads to quasi-stationary light field that undergoes a 
slow diffusive spread. The underlying physics of this process is identified as pulse 
matching of probe and control fields. We then show that spatially modulated control- 
field amplitudes allow us to coherently manipulate and compress the spatial shape of 
the stationary light pulse. These techniques can provide valuable tools for quantum 
nonlinear optics and quantum information processing. 

Key words: electromagnetically induced transparency, slow light, quantum 
information processing 



1 Introduction 



A very promising avenue towards scalable quantum information systems is 
based on photons as information carrier and atomic ensembles as storage 
and processing units 0. While a number of techniques for reliable transfer 
of quantum information between light and atomic ensembles have been pro- 




)osed^_^^_^ 
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and in part experimentally realized over the last couple of years 
[lol |. the implementation of quantum information processing in 



^ We dedicate this paper to Bruce W. Shore, one of the fathers of the theory of 
coherent processes in atomic systems, at the occasion of his 70th birthday. 
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these systems remains a challenge. This is because deterministic logic opera- 
tions require strong nonlinear couplings between few photons or collective ex- 
citations corresponding to stored photons. To achieve these strong couplings, 
long interaction times and tight spatial confinement of the excitations are 
needed. Even if long-range interactions between stored photonic qubits are 
employed as e.g. in the dipole-blockade scheme 0, 3 ^ tight confinement 
is needed to reach sufficiently high fidelities. 

In the present paper we discuss a method that could allow to manipulate 
the spatial shape of a collective excitation corresponding to a stored light 
pulse. It is an extension of the recently demonstrated technique to generate 



quasi- stationary pulses of light [IJ, ll5[ in electromagnetically induced trans 



parency (EIT) jl6| using counter-propagating control fields. In a light 
pulse was first stored in a delocalized state of an atomic ensemble by creating 
and adiabatically rotating the collective atom-light excitation, the so-called 
dark-state polariton from a freely propagating electromagnetic pulse into 
a stationary Zeeman excitation. The adiabatic rotation, which is accompanied 
by a decrease of the group velocity, is facilitated by reducing the intensity of 
the EIT control field. In the form of a pure Zeeman or spin coherence the 
excitation is stored and well protected from the environment for rather long 
times. At the same time it is however also immobile, thus preventing any 
manipulation of its spatial shape. In addition, the absence of any photonic 
component to the excitation prevents the use of nonlinear optical interactions 
for making such excitations interact. Regenerating a small photonic compo- 
nent of the polariton by means of a weak stationary retrieval field created by 
two counter-propagating lasers, a quasi-stationary pulse of light was created 
in the experiment of |l4l| in a second step. Through a mechanism known as 
pulse-matching Il8l | , the intensity of the regenerated stationary pulse fol- 
lows the oscillatory profile of the retrieve laser intensity. This process allows 
one to create a stationary excitation with a finite photonic component, i.e., an 
excitation with stationary, localized electromagnetic energy. One remarkable 
consequence of this effect is the possibility to enhance nonlinear optical pro- 
cesses jlH. Another interesting aspect of this technique is that despite the fact 
that the photonic component is at all times very small, the dark-state polari- 
ton becomes sufficiently mobile to follow the profile of the retrieval field. This 
provides a potential mechanisms to manipulate and control the spatial shape 
of a polariton while keeping most of its probability weight in well-protected 
spin coherences. 

In the present paper a detailed one-dimensional model of the generation of 
quasi-stationary pulses of light by counter-propagating lasers will be presented 
and its predictions compared to numerical simulations. It will be shown that in 
the weak-probe limit the dynamics of the regenerated light pulse is described 
by a set of coupled normal-mode equations 3] from which exact expressions 



for the temporal behavior of the pulse width can be obtained. It is shown that 
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for sufficiently large values of the optical depth (OD), these equations reduce 
to a simple diffusion equation, with the diffusion coefficient proportional to 
the group velocity divided by the optical depth. It is shown furthermore that 
control fields with spatially varying intensity profiles allow to manipulate the 
spatial shape of the stationary light pulse. Using a frequency comb of retrieve 
fields, a very narrow stationary mode proffie can be generated by a ffitering 
process. Alternatively using retrieve fields with an intensity difference that 
varies linearly in space in the region of interest will lead to a stationary field 
with a Gaussian spatial profile and an amplitude exponentially decaying in 
time. A linear dependence of the intensity ratio can be obtained e.g. by using 
paraxial retrieval beams with spatially displaced foci. Finally we demonstrate 
by numerical examples that moving the laser foci allows to shift and to com- 
press the stationary pulse of light. Although a quantitative estimate of the 
fidelity of such a compression process is not given here, this shows that sta- 
tionary pulses of light have a great potential for the manipulation of the spatial 
shape of stored photons. 



2 Stored-light retrieval with counter-propagating control fields 



2.1 Model 



Let us consider an ensemble of A-type three-level atoms with one excited level 
I a) and two lower levels \b) and |c). As shown in fig. la the transition |c) — |a) of 
the atoms is coupled to an external drive field of frequency Uc and wavenumber 
kc characterized by the Rabi-frequency fl+{z, t) = fl+Q{z, t) exp{i{kcZ — ujct)}. 
Respectively the transition |6) — |a) is coupled to a weak probe field of center 
frequency u and wavenumber k described by E^{z,t) = E^{z,t) exp{i{kz — 
ut)}. The excited state \a) decays radiatively with rate 7. All other decay 
and dephasing processes are neglected. Under conditions of two-photon reso- 
nance, i.e. for 6 = Ucb — — Uc) = 0, Ucb being the resonance frequency of 
the lower-level transition, the control field creates electromagnetically induced 
transparency (EIT) for the probe field. Associated with this is a reduction 
of the group velocity f gr of a light pulse within a certain frequency range close 
to two-photon resonance: 

Ugr = c cos^ 9, tan^ 9 = ^ . (1) 

+ 



Here g = d/h^hwahl'2(-o is the probe-field coupling strength proportional to 
the dipole matrix element d of the \a) — \h) transition, and is the density 
of atoms. An adiabatic rotation of ^ = 9{t) from a value close to zero to 7r/2 
leads to a slow-down and eventually to a full stop of the probe pulse, which 
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is associated with a transfer of its quantum state to a delocalized collective 
spin (Zeeman) excitation. An adiabatic rotation of the mixing angle from 7r/2 
back to zero (or some other value different from 7r/2) at a later time leads 
to the retrieval of the stored light pulse 0, H, |^ which then propagates in 
the original direction with a group velocity determined by eq.(l). If for the 
retrieval a drive beam is used with a different direction, the stored light is 
emitted into a direction determined by phase-matching |2ol . 2ll |. 



(a) 



lb>- 




(c) 




|si+ + o_p 



Fig. 1. Level scheme for storage of light pulse by intensity reduction of control 
field 0+ in an EIT system (a) , and subsequent generation of a stationary light pulse 
with two counter-propagating components E± by applying two counter-propagating 
control fields (b). The total field distribution of the retrieved light pulse and the 
control fields is indicated in (c). 

A very intriguing variant of the retrieval process was suggested and exper- 



imentally demonstrated in Rather than using a single coupling field, 
two counter-propagating retrieval beams of the same frequency and intensity 
were used. As indicated in fig. lb this leads to the generation of two counter- 
propagating probe fields Fj±{z, t) which form a quasi-stationary standing wave 
pattern as indicated in fig.lc. As shown in 0,0] the intensity profile of the 
retrieved field \E+ + 
that of the control field 
identical to the stored light field. Due to a process known as pulse matching 
17, 3 13' probe- field envelope tends to approach that of the control fields 



^ shows an interference pattern which is similar to 
+ The initial envelope of this pattern is 



with increasing time. As a consequence there is a diffusion-like behavior of the 
retrieved field envelope. In the following we want to theoretically analyze the 
underlying physics of this phenomenon within a one- dimensional model. This 
will then allow us to discuss a number of interesting generalizations. For this 
we consider the interaction Hamiltonian in a slowly-varying time frame 



H = -h dz 



Aacc{z,t) + {5 + A)aaa{z,t) + (Q{z,t)aac{z,t) + h.a.], 



+g E{z, t)aab{z, t) + h.a. 



(2) 
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where A = ujac — and b = Ucb — (cu — Uc) are the one- and two-photon 
detuning. The a^jy's are continuous versions of the slowly- varying (in time) 
single-atom flip operators 

(^^.Az,t) = Jim^^ ^M^W' = (3) 

In the above sum, the atom index i runs over all atoms with positions within 
the interval {z—Az/2, z+Az/2). Q{z, t) and E{z, t) are the (positive- frequency) 
complex Rabi-frequency of the drive field and the dimensionless slowly- varying 
complex amplitude of the probe field respectively. Both can be decomposed 
into two counter-propagating contributions 



n{z, t) = n+o{z, t) e"'^' + n_oiz, t) (4) 
E{z, t) = E+{z, t) e'''^' + E_{z, t) q-'^^' . (5) 

Note that a fast oscillating term with the retrieve wavenumber kc was split off 
also from the probe field. Making use of the commutation relations 

aap{z,t),af,y{z' ,t) 

and considering the perturbative regime of a weak probe field, where we can 
set abi,{z, t) ^ N ^ const., we find the following Langevin equations of motion 
for the atomic operators: 

CTba = -{^{S + A) + j)a,a + igN {E+e''^' + E.e"'^^^) 

+2 (fi+oe^'=^ + n.oe-^'^^') ahc, (7) 

CTbc = -^b^Tbc + ^ (fi^oe"*'^^ + f^loe*'^') ^ba. (8) 

We have dropped the Langevin noise operator in the equation for the optical 
coherence aba associated with the decay rate 7, since we want to work in 
the adiabatic limit in which this term is negligible jl^]. The above equations 
suggest the decomposition of the optical coherence in two counter-propagating 
components aba{z) = al^\z,t)e^''''^ + al^~\z,t)e~^'^''^ . Substituting this into (7) 
and (8) and making a secular approximation, i.e. collecting terms with the 
similar oscillatory terms (e='=*'^'=^) and neglecting fast oscillating contributions 
(^Q±z2kcz^^ yields 



5{Z - z') (Sp^ aau{z, t) - 5au (^^.^{z, t) 



(6) 



4? = -{K5 + A) + lY^'^ + ^9NE± + ifiioCTfe., (9) 

abc = -iSabc + iQ^o^ia^ + ^^^o4a ^ • (10) 

In the following we assume fi^Q = Q±o- Note, that it is also possible to not 
make the secular approximation, as is discussed briefly in the appendix. 
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If the temporal changes of the slowly varying field amplitudes are slow com- 
pared to 7~^, we can adiabatically eliminate the optical coherence. Under 
these conditions we find 

^''^^ ~ z(A + 5)+7 ' ^ ^ 

which leads to the effective equation of motion for the spin coherence 



\ z(A + d)+7/ z(A + d)+7 

where Qq = Q'^q + l^^g- Eqs.(ll,12) are the main equations for the tempo- 
ral evolution of the atomic system. They describe the dynamics of the spin 
coherence adiabatically followed by the optical coherences. The second set of 
equations needed for the description of the system are the wave equations for 
the two probe field components E±. In slowly- varying envelope approximation 
and within the one- dimensional model considered here, they read 



dt ± cd^j E±{z, t) = -zAu E±{z, t) + tgai^l (13) 

Here Au; = uj — Uc, and free-space dispersion k = ui/c and kc = uJc/c was 
assumed. 



2.2 effective field equations in the adiabatic limit 



In order to solve the shortened wave equations for the probe fields coupled to 
the atomic spins, we now apply an adiabatic perturbation to the time evolu- 
tion of the atomic spin coherence, eq.(12). In lowest order of this expansion we 
ignore the time derivative of the spin coherence ate altogether. This approx- 
imation is however not sufficient, since it does not capture important effects 
such as the group velocity and is only valid if the characteristic pulse times 
are large compared to the travel time through the medium. In order to de- 
scribe finite group velocities, the first order correction needs to be taken into 
account. This yields 



gN(E+n+o + EM^f 



ni + isr 

5) +7 



gNin+odtE+ + n^odtE^ 



(14) 
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where we have disregarded time derivatives of f2±o and have introduced the 
complex parameter F = 'y+i{A+6). Substituting this result into the expression 
for the optical coherence aba and subsequently into the wave equations (13) 
eventually leads to the coupled field equations 



± + ^ 



(15) 



These equations can be written in a more transparent form by introducing the 
mixing angles 6 and (j) 

tan 9 = —2-, tan = — — -, (16) 
iLq |"+o| 

and by assuming a small two-photon detuning, i.e. ^ |5r|: 



dt + c cos' 98, j E+ = -i [Acu cos' 9 + S sin' 9) E+ 

+ sin^ 9 sin (^sin - cos (pdtE^ (17) 

— cos^ Y (f^E^ — cos (j)EJ^ , 

and 

(dt - c cos^ 9d,^ E^ = -i (Aw cos^ 9 + 5 sin^ 9^ E_ 

— sin^ 9 cos (j) ^sin (pdtE^ — cos (pdtE^ (18) 

+ cos^ Y '^(^^^ 'PE+ — cos (pE^ . 

One recognizes that the two field components propagate with an effective 
group velocity Wgr = c cos^ ^ similar to eq.(l). The first bracket on the right 
hand side of eqs.(17) and (18) represents a phase-mismatch, which however 
vanishes for a two-photon detuning chosen such that 

S = -Aucoe9 ^ -Au^. (19) 

c 

If |fgr| <^ c the two-photon detuning is very small and does not lead to a 
violation of the EIT condition. In the retrieval process the probe field will 
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build up with a center frequency such that the phase-matching condition (19) 
is fulfilled. 



2.3 normal modes and pulse matching 



The structure of the two field equations (17) and (18) suggests the introduction 
of the two normal modes 



Es = cos (j) + siYKj) E_, Ed = sincj) E+ — cos (j) E_, (20) 

which we denote as sum and difference normal mode. In terms of these modes 
the propagation equations read 



dt + fgr cos{2(f))dzj Es + fgr sm{2(j))dz Ed = 

+Vgr{d,4>) {sm{2<j))Es - cos{2<P)Ed), (21) 

dt - ccos(2</.)9,j Ed + csm{2(P)d,Es = -^Ed 

-c[d,(j)) (cos(2(/))^5 + sm{2(j))ED) ■ (22) 

Here we have assumed that the mixing angles 6 and (p can be space dependent 
but are constant in time. At the same time, in lowest order of the adiabatic 
expansion, the atomic spin coherence follows the evolution of the sum normal 
mode Es, hence we find from eq.(14) 



abc = -VNtanOEs. (23) 

One recognizes from (21) and (22) that apart from the coupling between the 
normal modes Es and Ed, the difference mode Ed is strongly absorbed due to 
the term g'^N/T on the right hand side. As a consequence the amplitudes of the 
retrieved fields approach a configuration where Ed 0, i.e. a configuration 
where the probe amplitudes match those of the drive fields: 

|±^cot,,^|l±. (24) 
This phenomenon called pulse-matching is well-known for EIT systems [1 

SHI. 



One recognizes that for almost identical control fields, i. e. f2+ ~ the group 
velocity of the sum mode Es and difference mode Ed are very small due to 
the cos(20)-term. It is zero if the strength of the two counter-propagating 
control beams is exactly the same. If ^2+ is bigger than the sum mode 
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moves to the +z direction and vice versa. Finally even if the ratio of the 
control field envelopes is spatially constant, i.e. if dzcf) = 0, a small amplitude 
of the difference mode will be generated out of the sum mode due to the term 
proportional to dzEs until Es is constant in space. This coupling will give rise 
to the slow spatio-temporal evolution discussed in the following sections. 



3 Quasi-stationary pulses of light from spatially homogeneous re- 
trieval beams 



Let us first consider the case of two spatially homogeneous control fields with 
equal intensity, i.e. 

cos(20) = 0, dz<p = 0. (25) 

This can be realized e.g. by two laser beams of the same intensity with a 
negligible curvature of the phase fronts, i. e. in the plain wave regime. In this 
case the propagation equations for the sum and difference mode (21) and (22) 
simplify to 

dtEs = -v^AED, (26) 
dtED = -cdzEs-^ED. (27) 



3. 1 Adiabatic elimination of difference normal mode and diffusion equation 
for resonant probe fields 



Let us consider the case where the drive field detuning A is chosen such that 
the probe fields are resonant, i.e. A + 6 = 0. Since for an optically dense 
medium the phase matching condition (19) requires only a very small two- 
photon detuning, this is essentially equivalent to the case of resonant drive 
fields. Then F = '~f + i{A + 6) = 7 and eq.(27) shows that the difference normal 
mode is damped with a rate g'^N/^y = c//abs, where /abs is the absorption 
length of the medium in the absence of FIT. For a medium with sufficiently 
high optical density OD = I /labs the absorption length is typically on the mm 
scale and thus the decay time is on the order of a few picoseconds. The typical 
pulse times in light storage experiments are however much larger. Thus an 
adiabatic elimination of the difference normal mode, i.e. neglecting dtED as 
compared to {g'^N/'j)E£i, seems justified. Such an ehmination leads to 

d 

Ed = —labs^Es- (28) 
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Making use of this approximation we find for the sum normal mode Es a 
simple diffusion equation 



dt ^ 



abs 



Es 



92 



(29) 



In the retrieval process two counter-propagating probe field components are 
created with an initial envelope given by the stored spin excitation. These 
components then undergo a diffusion process with a diffusion constant D = 
"WgrLbs given by the product of group velocity and absorption length. This is 
illustrated in fig. 2 where false-color images show the two field distributions 
E+{z,t) and E_{z,t) for a storage process followed by a partial retrieve with 
two homogeneous, counter-propagating drive beams of equal amplitude. The 
data are obtained from a numerical solution of the wave equation (13) as well 
as the full set of atomic density matrix equations in secular approximation. 
The predicted diffusive behavior is nicely reproduced. 
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Fig. 2. Numerical simulation of storage and partial retrieval of a light pulse, where in 
the retrieval process two counter-propagating drive beams with equal and spatially 
homogeneous intensities have been applied. For the storage and retrieval we have 
used cos2 d+ = 0.5 * [1 - tanh(0.1 * {t - 65))] + 0.5 * 1/3 * [1 + tanh(0.1 * {t - 300))], 
where tan^ 9^. = g'^N/Q'^^. For cos^ 0„ only the second term was used. The color 
code represents the amplitude of the forward (£'+) and backward {E-) propagating 
components of the probe field. The parameters used for the numerical simulation 
are: 6 = A = Aui = and 7 = 1 and the initial width of the Gaussian wave-packet 
was Az{t = 0) = 10. 
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In the diffusion process the width of the probe field distribution as well as 
that of the collective spin excitation (see eq.(23)) increase according to 



Az'^it) = Az^to) + 2D{t~to). (30) 

Associated with this is a decrease of the excitation density. Since in a diffusion 
process the spatial integral of the field is constant but not the integral of the 
square of the field, representing the number of photons, there is also a (non- 
exponential) decay of the total number of excitations. After the control fields 
are switched on again, the sum mode has a Gaussian shape with width Az{0), 
and the total excitation, i.e. in the retrieved fields and the collective spin, 
evolves according to 

nUt) = ntot(O) (31) 
Az^{0) + 2Dt 

Thus in order to have negligible losses, the time over which a stationary pulse 
can be maintained is hmited by 

*«^^^. (32) 

U ^gr^abs 

which is exactly the characteristic time for the spread of the initial wave- 
packet. 



3.2 Small optical depth 



If the optical depth of the medium is small, the adiabatic elimination of En 
may be no longer justified. It is still possible however to find analytic results 
for the moments of the stationary light field. First of all one finds from eqs.(26) 
and (27) that like in the diffusion limit discussed above, the integral of Es is 
a constant of motion since E£){z = ±oo) = 0: 

Assuming an initially symmetric spin excitation around z = 0, the width of 
the retrieved light beam is given by the second moment of the sum mode 
d{t) = Az'^it) = dz Es{z, t)/A, which is coupled to the first moment of 
the difference mode g{t) = dz z ED{z,t)/A: 



-d{t) = 2v^,g{t), 



dt 



git) 



''abs 



-g{t) + c. 



(34) 
(35) 
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The solution of these equations can easily be found and reads 

dit) = diO) + 2Dt + (l - ^] (e-^*/'-= - l) . (36) 

C y tabs / 

One recognizes that a small absorption length /abs only affects the short-time 
evolution. In fig. 3 a comparison between the analytical prediction for 6{t) 
obtained from eq.(36) with Vgr — > fgr(^) and a numerical simulation of the 
Maxwell-Bloch equations is shown. Apart from a short initial time period, 
where due to the time dependence of fgr(i) non-adiabatic couplings lead to 
small deviations, there is a nearly perfect agreement between analytic predic- 
tion and numerical simulation. 



110 




105 1 ( ( , , , 1 

310 312 314 316 318 320 322 

t[1/(g N"2)] 

Fig. 3. Comparison between analytic prediction for Az^(t) of the stationary 
wave-packet of light from eq.(36) with a numerical simulation of the coupled 
Maxwell-Bloch equations. The data are based on the simulation of fig. 2. The in- 
sert shows the behavior on a larger time scale. 

3.3 Non-equal but constant drive intensities 

If the intensities of the counter-propagating retrieve beams are not equal but 
constant in time and space, one has 

cos(20) ^ 0, 9^0 = 0. (37) 

As a consequence the equation of motion of the sum normal mode attains, 
after adiabatic elimination of the difference mode, a finite drift term 

(dt + cos(20)9,) Es = %iabs sin' (20) d^Es. (38) 

Transforming into a moving frame with z' = z — Vg,- cos(20)t and t' = t 
leads again to a diffusion equation with a modified diffusion constant D = 
Wgr^abs sin'(20). Thus in the case of non equal, but constant drive intensities. 
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the diffusive behavior of the quasi-stationary hght is superimposed by a drift 
motion with a small velocity WgrCos(20). This can be understood in a very 
intuitive way. If one of the drive fields is stronger than the other one, Raman 
scattering occurs with higher probability into the probe mode co-propagating 
with the stronger drive field which causes a drift motion of the quasi-stationary 
wave-packet. 



4 Stationary pulses of light generated by spatially modulated re- 
trieve fields 



In this section, we discuss two techniques to manipulate the shape of stationary 
pulses with the ultimate goal of confining the pulse to very short spatial dimen- 
sions. Initially the pulse is stored as a spin coherence with a spatial envelope 
that extends over many wavelengths. We have seen in the previous sections 
that the (partial) retrieval of the stored light pulse by counter-propagating 
control fields leads to quasi- stationary pulses of light. The shape of these sta- 
tionary pulses is determined by the envelope of the initial spin coherence as 
well as the control-field envelopes through the mechanism of pulse-matching. 

This suggests two different mechanisms to manipulate the shape of the sta- 
tionary pulse of light. In the first method the atoms are illuminated by a 
frequency comb, i.e. with control beams that have multiple frequency compo- 
nents of equal intensity in the forward and backward direction. In this way 
many corresponding frequency components are generated for the signal field. 
These components interfere to create a very sharp spatial envelope, which is 
matched to the sharp spatial envelope of the control field, potentially confined 
over only a few wavelengths (see Fig. 4). 

A second method employs a spatial modulation in the difference of the for- 
ward and backward retrieve intensities. In the last section we have seen that 
unequal retrieve intensities can lead to a drift motion of the stationary field 
with an effective group velocity Vgj. cos{2(j)). If e.g. Wgi.cos(20) would be nega- 
tive for positive values of z and positive for negative values of z, the associated 
drift motion would tend to spatially compress the stationary field. As will be 
shown this can compensate the diffusive spread found in the last section. This 
situation, indicated in fig. 5, can be achieved when (and in general also 6) 
are made z dependent. 
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Fig. 4. (a) Atomic level structure with multiple-frequency-component 
control field U{z,t) = ^^0^6"*^'='^*"^/'^^ and generated signal field 
E{z,t) = Yl,k-^k[z^t)e~'^^''^^~^/^\ (b) Generated signal field envelope, show- 
ing tight localization due to constructively interfering frequency components 
(frequency comb). 




v cos(2(^) 



Fig. 5. Paraxial retrieve lasers with spatially varying, Gaussian intensity profiles 
determined by the position of the focal points z+ ^ Z- create a spatially varying 
effective group velocity. 

4-1 Shaping of stationary light pulses using the optical comb technique 

Let us first consider the case when the atomic sample is irradiated with sev- 
eral counter-propagating control fields, with detunings A^, and (complex) 
Rabi frequencies fl±k, such that the slowly- varying total Rabi frequency reads 
^{z,t) = Ek (fi+fce-^^^(*+"/'^) + (]_fce-^^^(*-"/'^)). Afc = tOk-uJab is the detuning 
from the atomic resonance. For simplicity we here consider a degenerate level 
scheme, i.e. Uab = ojac- We assume equal intensities of the corresponding for- 
ward and backward components = |^-/c|- Corresponding to these driving 
fields are signal field slowly varying amphtudes E±k{z, t), so that the total sig- 
nal field envelope is E{z, t) = Ek {E+k{z, t)e-^^fc(*-^/^) + E_k{z, t)e-^^fe(*+^/'=)). 

Assuming IA^^qI ^ I^/ct^oI so that we can ignore coupling of the frequency 
components through off-resonant processes, we may also expand the optical 
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polarization as aba{z,t) = Efc ((tJ+''^(z, t)e-*^^(*-^/'^) + t)e-^^^(*+^/=)) . 

Within the weak-probe and secular approximation and assuming two-photon 
resonance of all associated pairs of fields we obtain the following equations of 
motion 



ajf ) = - (z Afe + 7) ajf ^ + igNE±, + iQ^k^Tbc, (39) 

= + E l^^xAt'^ + ^^-A:'^] ■ (40) 

k 

Here again the Langevin noise operators associated with the spontaneous de- 
cay from the excited state in eqs.(39) have been neglected as they do not con- 
tribute in the adiabatic limit. The atomic polarizations 0"^^'^'' drive the probe 
field components E±k through the shortened one-dimensional wave equations 

(ifc) 



dt±cd,)E±k{z,t)=tgal,r- (41) 



Solving eq.(39) adiabatically for o"^^'^'' yields 



(±fc) WN iQ±k 



where we have introduced the notation Ffc = 7 + i^k- Substituting this result 
into eq.(40) for the ground-state coherence leads to 



^bc = -2^ f C^bc-gJy}^ :^ • (43) 

k ^ k 

Letting Aq = and |Afc^o| ^ 7? I^ifcl, we find that the spin coherence is 
driven only by the resonant fields, i.e. 

d,. - -%bc - ,Ar^lk^±!I±i^k^, (44) 

7 7 

where Qq = |fi+oP + |f^-oP- In a similar way as in Sec. 2. 2 we can solve this 
equation in first order of an adiabatic expansion. This yields 

ab. = -^iV ^+°^+o + ^^ol^-o ^ (ni,d,E^o + n^odtE^o) ■ (45) 

Substituting this result first into the equations for the resonant fields E±q 
leads to similar equations as in Sec. 2. 
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'dt ± ccos^ ed^E^^ = ±^fi;g(f]_o9iE+o - ^+^dtE.^) (46) 





Taking into account that = these equations can be written in a 

simpler form introducing sum and difference normal modes Es = {Q\qE+o + 
^%E^o)/no, Ed = (fi-o-E'+o — i^+oE^o) /Qq: 



dtEs = -v^AED, (47) 

dtEo = -cd,Es-^En. (48) 

7 

The second term on the right hand side of eq.(48) leads to a fast decay of 
the difference mode, such that in the long-time limit the resonant probe-field 
amphtudes are matched to the corresponding drive-field amplitudes 

Q.oE+o = ^+oE-o = (49) 

gN 

As in section 2, the sum normal mode undergoes a diffusion process under 
conditions when the difference normal mode can be adiabatically eliminated. 
As a consequence the spin excitation, which adiabatically follows the sum 
normal mode, does the same, i.e. 

dtabc{z,t) = Ddlabciz,t), (50) 

with D = fgr/abs- 

For the non-resonant probe-field components -Eifc, k ^ one finds the short- 
ened wave-equations 

(dt ± cd^ E^,{z, t) = -^E^, - ^a^c. (51) 

V / i A: J- A; 

The second term on the right hand side does not depend on E±^.. As a conse- 
quence, assuming a sufficiently dense medium, the off-resonant probe ampli- 
tude may be adiabatically eliminated leading to 

E±k = ~'^^bc- (52) 

Noting that a similar relation holds for the resonant components in the long- 
time limit, eq.(49), we finally arrive at 
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gN 



- E (f^+.e-^^(*-/^) + f]_,e-^^(*+^/^)) (53) 



o-fec(^,i) 



-fi(z,t) 



gN 



Thus the electric field envelope E{z,t) becomes matched to the control field 
envelope Q{z, t) modified by the spin coherence abc{z, t). This allows to control 
the stationary pulse shape through control of the retrieve amplitudes and 
phases. 

Note that the condition for this analysis to hold is kprn/Lhs ^ niax{Afc}/7. 
Taking the length of the initial spin excitation 4pin to be such that the pulse 
just fits inside the medium /spin ~ this condition implies that the maximum 
frequency detuning A^ax for which the stationary pulse can adiabatically fol- 
low the control field through pulse matching, is given by A^ax ~ 7 OD, where 
OD = = l/Uhs is the on-resonance optical depth. Thus spatial features 
as small as the optical depth can be imposed on the stationary pulse through 
the frequency comb technique. 

It should be noted, however, that the generation of spatially narrow stationary 
fields by means of the frequency-comb technique is a filtering process rather 
than a compression of excitation. In fact the total number of probe photons 
created by a frequency comb is much less than in the case when only the 
resonant components Q±.q of the retrieve laser are present. The excitation 
density at the center of the stationary photon wavepacket is the same in both 
cases, while in the wings it is substantially smaller for the case of the frequency 
comb as compared to the case of homogeneous retrieve beams. 



4-2 Shaping of stationary light pulses using a spatially varying group velocity 



Let us now discuss the second method indicated in fig. 5 in detail. Assuming 
again single and two-photon resonance and an optically thick medium, we can 
adiabatically eliminate the difference normal mode from (22). This yields 

Ed-- sm{2<P)ky,AEs - cos(20)/abs (5,0) Es, (54) 

where we have assumed that changes only little over the absorption length 
Zabs and thus \{g'^N/'y)\ ^ 1(9,0) sin(20)|. Substituting this into (21) gives 









AiEs 




DEs 



(55) 
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where D = WgrLbs is the diffusion constant introduced before, and the coeffi- 
cients Aq and Ai read 



Ao- 



Ai 



-V, 



gr 



—V. 



{d,(t)) sin(2</.) + 2/abs(9,</.)' sin2(2(/.) 
-/abs(5.0)'cos2(20) - {diet)) cos(20) sin(20) 

l+4/abs('9.0) 



(56) 
(57) 



The constant term proportional to Aq in eq.(55) can be removed by the sub- 
stitution 

Es = Es exp{Aot}, (58) 
which results into a Fokker-Planck equation for Es'- 



d_ 

di' 



E, 



d_ 

dz 









AiEs 


dz'^ 


DEs 



(59) 



In the following we will discuss the spatio-temporal evolution of Es resulting 
from this equation. 

Non equal drive fields lead to an effective group velocity Vg^ cos (20) for the sum 
normal mode. If this group velocity is tailored in such a way that it is negative 
for positive values of z and positive for negative values of z, there is an effective 
drift towards the origin. This force may compensate the dispersion due to the 
absorption of large- /c components of the probe field found in section 3. We 
thus consider as the simplest example the special case of a linearly varying 
intensity difference of the two drive fields with a constant sum = const.: 



cos(20) 



sin (20) ^ 1 . 



(60) 



This situation is realized e.g. if the two control fields are paraxial, Gaussian 
laser beams with focal points at z± ^ ±21. The linear approximation is of 
course only valid for \z\ <C /. In this case eg. (59) turns into the Fokker-Planck 
equation of the Ornstein-Uhlenbeck process [2J| for which exact analytic so- 
lutions are known 



d_ 

dt 



Es 



d_ 

dz 



z ~ 






VgrjEs 


dz'^ 


DEs 



(61) 



Here we have neglected contributions proportional to /abs/^ as compared to 
unity. The Ornstein-Uhlenbeck process has a stationary Gaussian solution 



with width \^TI^. Noting that now Aq 
eq.(58) in the long-time limit: 



Es{z,t) = Es{z) exp|Aot} — > exp 



2U 



-gr 



exp 



/2l this gives with 



abs 



(62) 
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Fig. 6. The prefactor of the drift term for the configuration of to separated Gaussian 
beams. The focal points of the two beams are located at z± = ib20^abs- The linear 
regime under consideration in this section is clearly visible. 

The use of retrieve lasers with non-equal and spatially varying intensities thus 
acts like an effective cavity for the probe field with a ring-down time l/vg,- 
given by the time a photon travels between the intensity maxima of the two 
drive lasers. 

The initial- value problem of the Ornstein-Uhlenbeck process can be solved by 
making use of the eigensolutions {^n{z), A„} of the corresponding backward 
equation 
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D I dz ''^ ' 



X 



D 



(63) 



Eq.(63) is the differential equation of the Hermite polynomials and thus 
the eigenvalues A„ and eigenf unctions $„(2;) read 



An = n—, 



$„,(^) = (2"n!)-i/2i7„ 



riG {0,1,2,...}, 



/ z 



\^/2l Zabs' 

The general solution of the initial value problem then reads 



(64) 
(65) 



^2 



V2"+^n!7r//abs [ 2//abs 



X 



The coefficients Cn are determined by the initial field £'5(^,0): 

c„ = rdzEs{z,0)HJ-^) {2-n\)-'/'. (67) 
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It is interesting to note that, apart from the additional overall damping term, 
eq.(66) is very similar to a damped harmonic oscillator with oscillator length 
a// /abs- If the stored light pulse is a Gaussian and if the separation of the foci 
of the two retrieve lasers is chosen such that the width of the stored light pulse 
is less than the effective oscillator length a// /abs, only the fundamental mode 
$0 gets excited in the retrieve process. In this spatially constant field 

distribution is created. The field has however a finite lifetime determined by 
the decay rate 7efr = fgr//. As a consequence the total excitation decays in 
time according to 

ntot{t) =ntot(0) exp|-t;grt//|. (68) 

In order to have negligible losses, the time over which the stationary pulse can 
be maintained is limited by exactly the same expression as in the diffusion 

case 

(69) 

Vgj- "^gr^abs 

If the separation of the focal points in the retrieve process is much smaller 
than /gpin/ /abs the generated stationary light pulse has a much narrower width 
than the original spin excitation (/spin)- In this case a large number of higher- 
order Gauss-Hermite modes is excited however (see eq.(67)), which decay much 
faster than the fundamental mode. Thus as in the case of the frequency comb 
very narrow spatial distributions of the stationary field can be created, however 
only through a filter process. 

In fig. 7 a numerical simulation of the the retrieval using two control beams 
with separated foci is shown. Initially a Gaussian probe pulse is stored in a 
collective spin excitation. The center of the stored probe pulse is in the middle 
between the two foci, indicated by the two white lines and the width of the 
pulse is on the order of a/T/^. Thus mainly the fundamental mode $o is 
excited, which can be seen from the constant spatial shape of the retrieved 
wave-packet. 



5 Spatial compression of stationary pulses of light 

We have seen in the last section that the use of retrieve lasers with spatially 
modulated intensities does allow the generation of stationary light pulses with 
very narrow spatial shapes. The underlying process is however a filter pro- 
cess and thus accompanied either by reduction of the photonic component in 
the polariton or by large losses. Nevertheless both techniques open interesting 
possibilities for the spatial compression of a stored photon with small losses. 
If a stationary light pulse is created e.g. with separated foci of the drive lasers 
as explained in the previous section, and the distance between the focal points 
is reduced as a function of time in an adiabatic way, the spatial width of the 
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-50 -30 -10 10 30 50 -50 -30 -10 10 30 50 

z[c/(gN^''2)] z[c/(gN^^')] 



Fig. 7. Storage of a Gaussian pulse and subsequent retrieval with two control 
beams with spatially changing intensity difference. The maxima of the intensities 
of Q± are indicated by the two white lines. Close to the midpoint between these 
lines the normalized intensity difference (Ifi+oP — |ri-oP)/^^o varies linearly with 
z. The generation of fields with constant spatial shape is apparent. The parame- 
ters are that of fig. (2). In addition Gaussian beam profiles 0,±{z) = l/w±{z) with 
w±{z) = \J\ — 2{z — z±)/tt have been assumed. 

stationary pulse follows. This results in a spatial compression of the probe 
excitation. As the effective decay rate of a stationary pulse 7eff increases with 
decreasing pulse width, the control fields should be switched off immediately 
after the compression. Fig. 8 shows a numerical example for such a process. 
After retrieval of a stored pulse with separated foci, whose position is indi- 
cated by white lines, the distance between the focal points is reduced. One can 
see very clearly that the field density as well as the density of the spin excita- 
tion is substantially increased in this process. This suggests that by changing 
the spatial profile of the control field in time, either using the optical comb 
technique or displaced foci, spatial compression can be achieved. 

In fig. 9 the temporal evolution of the peak density and the total excitation 
(i.e. photon number of the stationary field plus spin excitations in the atomic 
ensemble) are shown for the example of fig. 8. One recognizes that when the 
peak density increases the total number of excitations starts to decrease very 
rapidly. Clearly an optimization is needed to maximize fidelity and compres- 
sion. In order to find optimum conditions and to estimate the maximum 
possible compression from the theoretical model, it is necessary to include 
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-50-30-10010 30 50 -50-30-10010 30 50 -50-30-10010 30 50 

z[c/{gN^^^)] z[c/{gN^''2)] z[c/{gU"^)] 

Fig. 8. Retrieval of a stored pulse and spin excitation using drive fields with 
separated foci and subsequent reduction of their separation. The position of 
the focal points is indicated by white lines. One clearly recognizes a compres- 
sion of the probe field associated with an increase of the field density. The pa- 
rameters are the same as in fig. (7) except for 7 = 0.05. The foci move like: 
z±{t) = 20 =F 10 * 0.5 * [1 + tanh(0.0125 * {t - 700))] 



non-adiabatic couplings into the description. Finally, it is not clear if spatial 
compression of this kind can be used to enhance further nonlinear optical pro- 
cesses using the techniques of Ref. This analysis is beyond the scope of 
the present paper and will be discussed in detail elsewhere. 




Fig. 9. Peak excitation density (solid line) and total excitation (dashed line) as 
function of time for the example of fig. 8. All quantities are normalized to their 
value at t = 200{gVN)~'^. 



22 



6 Summary and outlook 



In the present paper we have discussed the generation and coherent control 
of stationary pulses of light by storage of a light pulse in a collective spin 
excitation via EIT and subsequent partial retrieval of this excitation with 
counter-propagating retrieve lasers. We have shown that for equal and spatially 
homogeneous intensities of the forward and backward retrieve beams a quasi- 
stationary wavepacket of the probe light is generated with an initial envelope 
given by the spin excitation. For an optically thick medium the dynamics of 
the wavepacket is described by a diffusion equation with a diffusion constant 
given by the product of group velocity and absorption length without EIT. The 
physical mechanisms for the diffusion of the stationary pulse is the well-known 
phenomenon of pulse matching of probe and retrieve field components in EIT. 
We have shown furthermore that spatially modulated retrieve lasers can be 
used to manipulate the shape of the stationary light pulse and in particular 
to spatially compress the excitation. Making use of a frequency comb for the 
retrieve fields a very narrow spatial distribution of the probe field can be 
generated. Likewise the use of retrieve fields with spatially varying intensity 
difference can lead to a narrow non-dispersing but exponentially decaying field 
distribution. In both cases the narrow field distribution is however created by 
a filtering process. Thus these techniques can not straightforwardly be applied, 
e.g. to applications in quantum nonlinear optics. However, as demonstrated 
with some numerical examples, if the retrieve field distribution is initially 
matched to the stored spin excitation and its shape is modified in time in 
an adiabatic way, the stationary light pulse and thus the stored excitation 
can be spatially compressed. Although we have not analyzed the fidelity of 
the compression process quantitatively and have not optimized it, the present 
paper shows the potentials of coherent control of stationary light pulses for 
quantum nonlinear optics and quantum information processing with photons 
and atomic ensembles. 
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A Multi-component Spatial Coherence Grating 



In this appendix, we analyze the stationary pulse solutions from the point of 
view of spatial coherence gratings jisl. l26l|. and we show that multiple atomic 
momentum components can be taken into account. These arise due to multiple 
scattering of photons in the forward and backward directions, resulting in 
distinctly different atomic susceptibilities. For stationary atoms, such as in 
cold atomic samples, these multiple scattering momentum components can be 
populated and preserve their coherence. In contrast, for warm atomic vapors, 
the rapid random motion of the atoms and their collisions results in a very 
rapid decay of spatial coherences with period equal to or shorter than the 
optical wavelength. 

As discussed in section 2.1, associated with the forward/backward propa- 
gating fields are slowly varying optical coherences (polarization) with slowly 
varying envelopes P±{z,t), so that the total polarization can be written as 
\/Naba{z,t) = P+(z, t)e*'^=^ + P_{z,t)e~^^''^ . Similarly, the ground-state spin 
coherence is defined as S{z,t) = ^/Nabc{z,t). Letting the wave- vector mis- 
match be AK = kc-k, and defining S± = E±e^^^^'' and V± = P±e^*^■^^ 
the equations of motion for the fields can then be written as 

{dt + cd^)S+{z,t) = iAKc£+ + igVNV+ (A. la) 

{dt - cd,)8_{z,t) = lAKcS. + igVNV- (A.lb) 

while the atomic equations of motion are (setting A = 6 = 0) 

dtV+ = -jV+ + iQ+oS + igVN£+ (A.2a) 
dtV- = -7P_ + in^oS + igVNS. (A.2b) 
dtS = in\^V+ + in%V- . (A.2c) 



These equations show that the counter-propagating control fields f2±o in- 
duce a coupling between and mediated through the spin coherence 
S. This leads to the formation of new eigenmodes of propagation, where as 
we show below, there is one mode that is very rapidly decaying while the 
other decays very weakly in the large optical depth limit. This phenomenon 
is analogous to the "pulse matching" phenomenon, as first described by Har- 
ris ^3]- These equations can be used to obtain the susceptibilities, given by 
^Xaa' = g^/NVj£,., where a, a' = ±. 

We now contrast two approaches to compute the susceptibilities, one approach 
in which the secular approximation is made and one where it is not. Writing 
the polarization as V{z^t) = Y.n'P2n+i{z,t)e^^'^^~^^^''''^ , and the spin wave be 
<S{z,t) = Y.nS2n{z,t)e'^''^^''^ , we have the effective Hamiltonian 
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+ ^(n+o52nP2n+l + ^-0<S2n'P2n-l) + h.C. 



The equations of motion for the fields are given by 



(A.3) 



{dt + cd,)S+ = iAKc8+ + igVNVi (A.4) 

{dt - cd,)£^ = lAKcS^ + ig^/NV^i, (A.5) 

while the atomic equations of motion are 

dtV2n+l = -7^2n+l 

+ iQ+oS2n + ^^-0'52(n+l) + W'/N {6^,0^+1 + Sn,^lS^i) (A. 6) 

dtS2n = -l{^l^V2n+l + ^%V2n^l) (A.7) 



The susceptibilites Xcrcr'i!^) can be computed from these coupled equations. 
Truncating the decomposition of the polarization at V±{2n+i)i we can easily 
show that = reproduces our earlier result of section 2.1, whereas letting 
n — i> oo leads to a different limit (see Fig. A.l). We thus find that there is a 
clear difference between the multi-component case and the case when only the 
zero momentum component of the spin wave contributes. These two situations 
correspond e.g. to cold atomic samples vs. hot atomic gases, so that based on 
the previous considerations we expect that the stationary light pulse will be 
very different in cold atomic samples. 

The susceptibilities obtained in the limit n — > cxd can be shown 2^ to be the 



same as those found from a coupled mode approach |27|. Maxwell's equation 
in a spatially modulated medium in ID is (in the frequency domain) 



c^— E(^, v) + + x(^, y)]E{z, v) = 0. (A.8) 

Putting E{z, v) = i?+(2;, u;)e*'''=^ + £'-(2;, u;)e"*'''=^, where u = u — Uc, coupled 
mode equations can be obtained by letting the susceptibility be given by the 
usual EIT susceptibility [H 



2g'N t^Tp 

where the control field amplitude is space dependent Q{z) = Q.^.^e''^"^ + 
Q-oe~''^''^ , and where T = •j — iu, and Tq = 70 — iuj (70 is the ground state 
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coherence decay rate, which we ignore for simphcity in this work). Using the 
Fourier expansion x{z,uj) = Y.nXniyj)^-^'^^''^ and using coupled mode analysis 
[irj l. we can compute the susceptibilities Xa,cr'{'-^)- The result of this calculation 
is compared to the result of the multiple spatial component approach outlined 
previously, in Fig. A.l. 

Implicit in the derivation of the susceptibilities with a coupled mode approach, 
is the assumption of stationary atoms, so that all intermediate atomic mo- 
mentum states are taken into account, in particular high spatial frequency 
coherences contribute to this result. All spatial components of the spin and 
polarization waves are taken into account, even those with momentum equal 
to multiples of the optical wavevector K. For warm atomic samples, these 
coherences rapidly decohere due to the random motion of the atoms and their 
collisions. Hence, only for cold atoms, such as those in a Bose-Einstein conden- 
sate, these spatial coherences may play a role and lead to interesting differences 
with stationary light phenomenon observed in thermal atomic vapors. 

Finally, we argue that for storage and retrieval of excitations, these large spa- 
tial wavevector coherences are mostly irrelevant. As can be seen from A. 6 and 
A. 7, the signal field (in modes £±) couples to the polarization components V±i 
with Rabi frequency gy/N, which in the "slow" light limit is much larger than 
the control field Rabi frequencies fi±o- Therefore, starting from a stored pulse 
in the zero- momentum spin coherence So, most of the amplitude is coupled 
to the signal field modes S±, while very little amplitude "leaks" to the higher 
momentum coherences. Thus coupling to higher momentum coherences does 
not lead to decay of the spin coherence, as discussed in ji^. 
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